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Abstract 

The determination of the density functions for products of random elements from 
specified classes of matrices is a basic problem in random matrix theory and is 
fV| also of interest in theoretical physics. For connected simple Lie groups of 2 x 2 

matrices and conjugacy and spherical classes a complete solution is given here. The 
problem/solution can be re-stated in terms of the structure of certain Hecke algebras 
attached to groups of 2 x 2 matrices. 
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1 Introduction 

00 , 

Let (G, K) be a symmetric pair where G is a semi-simple Lie group with finite center. For 
G compact let H c = 7~t c (G) be the algebra under convolution generated by the invariant 
measures concentrated on conjugacy classes in G, and 7i s = H S (G, K) be the convolution 
algebra generated by the invariant measures on spherical classes O a = KaK C G. It 
is elementary that 7i c and 7i s are commutative algebras with unit. These algebras arise 
naturally in random matrix theory, and theoretical physics. In fact, the product structures 
. ^ \ for generators of these algebras are given by the density functions of products of random 
matrices chosen according to invariant measures on conjugacy or spherical classes; and 
in string theory generic D-branes are localized along the product of (twisted) conjugacy 
classes of the Lie group. We will not discuss specific physical interpretations in this paper 
(except for a remark in §5), and the reader is referred to [Q], and references thereof, for 
a detailed discussion of the issues of interest in string theory. 

This work may be regarded as the initial attempt at understanding the structure of 
the Hecke algebras 7i c and 7i s and considers only the special case of 2 x 2 matrices. The 
support of the density function for products of two generators of TC C (SU (n)) is determined 
in [AW] and is exhibited by a set of linear inequalities on the Lie algebra of a maximal 
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torus. The result in [AW] is essentially a reformulation of a theorem about singular 
connections on a holomorphic vector bundle on a Riemann surface with marked points 
[B], or equivalently a theorem of Mehta and Seshadri in algebraic geometry. Neither 
theorem is applicable for the computation of the density function for the product of two 
generators of the Hecke algebra. 

A complete structure theorem for 7i c {SU{2)) is given in Theorem 12.11 below. For 
spherical classes attached to groups of 2 x 2 matrices we consider the symmetric pairs 



(££7(2),£(t7(l) x [/(!)), (£L(2,M),£0(2)) and (£L(2, C), ££7(2)). Theorem O gives 



density functions for products of spherical classes. In the final section we present some 
numerical results in relation to products of conjugacy classes in £77(2). Each conju- 
gacy class, with the induced metric, is a copy of £ 2 equipped with a metric of constant 
curvature. It is noticed (perhaps surprisingly) that if £ 2 is dicretized according to the pre- 
scription of Thomson's problem (minimizing Coulomb potential or intuitively "the best 
equally spaced distribution") then convergence to the predicted measure is much slower 
than if the points were chosen randomly. 

In [JM] the question of whether a product of conjugacy classes C ai , . . . , C an contains the 
identity element and/or is the entire group ££7(2) is studied. By successive applications 
of Theorem 12.11 of this paper one can recover the results in [JM] and in fact give a more 
precise version of it. This subject has not been elaborated on here. The argument in [JM] 
depends on the results in [AW], but the proofs presented here are self-contained and of a 
more elementary nature. 

The author wishes to thank A. Katanforoush for the numerical results, R. Szabo, M. 
Sheikh- Jabbari, S. Shahshahani and especially Mehrdad Shahshahani for many stimulat- 
ing discussions. 

2 Statement of main results 

A conjugacy class Cg C SU (2) is uniquely determined by the eigenvalues e ±id of a matrix 
in Cg. Each conjugacy class is a homogeneous space, and for Cg this measure n$ is 
uniquely normalized to be 4% sin 2 6 in accordance with Weyl's integration formula and 
the normalization vol. (££7(2)) = Air 2 . 

Theorem 2.1 Let \i a and \Lp be the invariant measure on the conjugacy classes C a and 
Cf3 respectively (regarded as singular distribution on G). Then \i a * \ip is an absolutely 
continuous conjugation invariant measure on G relative to the Haar measure and its 
density is given by 




for a - (3 < 9 < a + (3, 
for - a - 13 < 9 < -a + (3, 
otherwise, 



(2.1) 
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where < (3 < a < ix. 

Let K = S(U(1) x U(l)) C SU(2), then (G,K) is a symmetric pair of compact type 
[HI]. Orbits of the action of K x K on G via 

a — > kak', (k, k') G K x K, a G G 

are called spherical classes and are denoted by O a . O a is a homogeneous space and 
it is a simple calculation that |on| is constant on O a and uniquely determines it. Here 

a — ( ° n 0,12 ) G SU{2) is any matrix in O a . Each spherical class O a carries an invariant 



^21 a 22 /f 

measure A a . The total mass of the measure A a of a spherical class O a is 2|an], and we set 
r(a) = \a n \. 

In the non-compact cases (SL(2,R),SO(2)) and {SL(2, C), SU(2)) the real Cartan 
subgroup A is 

A = {a t =( e l ° L ) \teR} 
\U e ay 

A + = {a t G A | t > 0} parametrizes spherical classes. The Haar measure for the Cartan 
(polar) decomposition G ~ is 

/ fWg = c f f f f{k 1 ak 2 )5{t)dtdk l dk 2 (2.2) 

where c is a suitable constant (see [H2]), 8(t) = sinh 2 t, and e = 1 or 2 according as 
G = SL(2,M) or SL(2,C). The volume (area) of the spherical class O at is 

vol.(C at ) = [47r 2 csinht] e . (2.3) 

Theorem 2.2 The product formula for spherical classes in the three cases of symmetric 
pairs attached to groups o/2x2 matrices are: 

(A)- Let X a and Xj, be two (singular) spherical measures concentrated on the spherical 
classes O a and Ob respectively. Then \ a * A& is absolutely continuous relative to the Haar 
measure on SU(2) and given by 

le^lau&nlu for ^ CQ _ Cl < u < y/co + d, 
V / C 2_( M 2_ C0 )2 ^ 2/ ^ 

0, otherwise, 
where Cq and C\ are symmetric functions of a and b and given by 

Co = r 2 (a)r 2 (6) + (l-r 2 (a))(l-r 2 (b)), 
ci = 2r(a)r(6) v / (l - r 2 (a))(l - r 2 (b)). 
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(B) - Let and A at2 be the (singular) invariant measures concentrated on the spherical 
classes O at and O at2 in SL(2,M). Then \ a *\ a is spherical and absolutely continuous 
relative to the Haar measure, and for a continuous spherical function f on SL(2, M) we 
have 

f Slllll T 

Aa* A 6 (/) = 4c 2 7r 2 sinh £i sinh £ 2 / fir) _ — = dr , 

Ji tl ,t 2 V c 2 ~ ( c i _ coshr) 2 

where I tl> t 2 — \p2 — Hi h + ti], Ci = coshti cosh£ 2; and c 2 = sinhtx sinht 2 . 

(C) - Let X ati and A at2 be the (singular) invariant measures concentrated on the spherical 
classes O at andO a±2 in SL{2,C). Then X ati * X at2 is spherical and absolutely continuous 
relative to the Haar measure, and for a continuous spherical function f on SL(2, C) we 
have 

A atl * A at2 (/) = 32c 2 7r 6 sinh t\ sinh t 2 / f(r) sinh rdr , 

where J tljta = [t 2 - *2 + h] . 

It may be of interest to normalize the measures on the spherical classes to probability 
measures and determine the corresponding empirical measure of products. For such a 
normalization the density functions determined in Theorem 12.21 become 

1 u g 1 sinhr sinhr (2 5) 

2^k ^Jc\- (n 2 ^' ' v/cl - (ci -coshr) 2 ' ' 2 sinh t x sinht 2 ' 

Remark 2.1 Let 7i c and TC S denote the completions of TC C and Tt s in the weak topology. 
Then in all cases considered 7i c and 7i s contain the corresponding L 1 space as a dense 
ideal. Furthermore, by the above analysis, for every pair of generators fj,± and /U 2 , we have 
III */i 2 G L 1 . 



3 Proof of Theorem ED 

First we show that fi a * is an LP function for p < 2. The Fourier expansion of the 
singular measure fi a is given by 

J2d P ^Wa)p(g)}, (3.1) 

peG* 

where G is the set of irreducible representations of G and the Fourier transform of a 
measure \i is defined as 

Pit*) = / p{9~ X )dp. 
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This series does not converge in the ordinary sense of convergence of series of functions 
since the measure p a is singular, however, it converges in the weak sense. The convolution 
product p a -k pp can be calculated from 

dpTr [p{p a )p(pp)p(g)] ■ (3.2) 

Since the character xp(d) is independent of g G C a , and is given by for P is the 

symmetric n th representation of G = SU(2), 

4 sin 2 a 4x P (C a )sin 2 a 

p{p a ) = Tr(p(^))/ = /. 

d p d p 

Applying the Plancherel theorem to p a * pp, we obtain 

ii ^ 1 12 16sin 2 asin 2 /3x P (C , a)Xp(^/3) , 2 

\\Pa*Pp\\ - 2^ ^4 d P 

= ^ — sin 2 (n + l)asin 2 (n+ 1)/?. 

This series converges absolutely and therefore p a * pp is a square integrable function. 
From the Cauchy-Schwartz inequality it follows that p a * pp is absolutely integrable. 

Therefore the measure p a * pp is a conjugation invariant function and it can be inter- 
preted as the (defective) density function for the space of solutions c to the equation 

abc = e, where a G C a , b G Cp. 

Since this density is conjugation invariant we represent it as the function v{6) on the 
maximal torus of diagonal matrices given by 

v{0) = ((p a *pp) *p e ){e). 

Let 1Z denote the left regular representation of G, and A a = An sin 2 a be the factor 
appearing in Weyl integration formula for conjugacy invariant functions. The Fourier 
transform of a function tjj on G at a representation p is defined as 

P(V') — / p(x~ 1 )ip(x)dx. 



Let dp denote the dimension of the representation p, Xp its character and G the space of 
(complex) irreducible representations of G. Taking Fourier transform and decomposing 
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the regular representation 1Z of G in the usual way we obtain 
u(0) = {{ix a * up)* no){e) 

= vo i fQ\ TrlZ ((^ * fip) * Me) 

= ^]7^y ^dpTip^* Hp) */^) 

peG 



_1 ^^A«^A,^ (a) ^ (/5) ^ WTr(/) 
peG p 



vol.(G) ^ rf p 

pSG 



Therefore 



,/,, • ' • ■■> < ■ '/) Xp{ a )xp{P)Xp{@) ,,,,,, 

'(f) = lD7rsiti asm pan S >, — j (A *) 

Up 

peG 

Irreducible representations of G = SU(2) are determined by dimension k > 2. The 
character of the dimensional representation pk is: 

' sin 6* 

i=o 

where e ±ld is the eigenvalue of the conjugacy class of g. Substituting in ()3.3j) we obtain 

„ sin(A; + l)a sm(k + l)3 sin(k + l)9 
v{6) = 16n sin a sin 8 sin 6 V — - - ^ — — ^- i — . (3.4) 

k=0 

Now let f(9) denote the function defined by (j2.1|) . Since / is an even function its Fourier 
expansion is of the form f(9) = ^ + J2 a n cosnx, and 

2 1 1 

a n = — T sin(n + l)asin(n + 1)3 sin(n — l)asin(n — 1)3\ 

7r n + 1 n — 1 

Substituting in the Fourier expansion we obtain 

>. .sin(m + l)a sinfm + 1)3 sin(m — l)a sinfm — 1)3. 

= 87rsinasin/3 > — - - - — - - - —)cosm6. 

/ — ' m + 1 m — 1 

m=0 

Using the elementary identity cos m^- cos(m + 2)# = 2sin(m+l)#sin6' the series becomes 
telescopic and simplifies to 

sin(m + 1)ol sin(m + 1)3 sin(m + 1)6 



f(9) = 167r sin a sin 3 sin 9 

m=0 

which is identical with ()3.4|) . ■ 



m + 1 
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Remark 3.1 In the above proof we made use of the fact that the measure p a * is a 
conjugation invariant function that can be interpreted as the (defective) density function 
for the space of solutions c to the equation 

abc = e, where a G C a , b G Cp 

and this density is given by ()3.3|) . For the case of finite groups this formula is well-known 
([S],p.68). 

Remark 3.2 It is possible to prove Theorem l2. ll without the use of harmonic analysis and 
by integral formulae similar to those used for the proof of Theorem 12 . 21 below. However, it 
appears that the above argument is possibly generalizable to products of conjugacy classes 
in compact connected semi-simple Lie groups, but the one based on integral formulae is 
not. 

4 Proof of Theorem 12.2 

Since the proofs of (B) and (C) are essentially the same computation we prove (A) and 
(C) only. Introduce coordinates on SU (2) by: 

(P, <P, V) — (_ J^Jj* V/I ^?r # ) ' (P> V> G P. 1] x [0. 2 ^ x P. ^ 

The Haar measure on SU{2) in the (p, tp, if})- coordinates is easily calculated by computing 
g~ 1 dg, a basis of left invariant 1-forms u>i,uj2,u}3 and then taking their wedge product to 
obtain: 

Q = uj\ A uj 2 A u 3 = 2pdpdtpdip. 

With this normalization vol(SU(2)) = 4-k 2 as before. 

Since both / and A a are K-bi-invariant, A a * f(x) is iT-bi-invariant and therefore to 
compute p a -k f{x) we can assume that x is of the form 

where t is a real number and w = se ia a complex number with t 2 + s 2 = 1. In (p, tp, ip)- 
coordinates on SU(2) we have 

p p2n pI'K 

\a*f(x)= f(yx~ 1 )dy = 2r(a) / / f(y(p,(p,ip)x~ 1 )dipdilj 

JOa JO JO 
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With x represented as in (|4.2j) we have 

w(-,-^.,ff/((_^^r)(4.r»^ 



Since / is spherical, it depends only on the norm of the (1, 1) entry of the above matrix. 
The square of the norm of the (1,1) entry is 

t 2 r 2 (a) + s 2 (l - r 2 (a)) + 2tsr(a)^Jl - r 2 (a) cos(^ + ip - a). 
Substituting t = r(x) and s = \/l — r 2 (x), the norm of the (1, 1) entry becomes 
\r(a)te itp + by/l - p 2 e~^\ = ^c + c x cos(y? + ^ -a). 

Therefore 

X a -kf(x) = 2r(a) \ \ f(y/ cq + Ci cos(ip + tp — a)) d(p dtp. 
The change of variable 



o jo 



(u,v) = (a/co + ~c\ cos((^ + if) - aj,ip), 
is a 2 to 1 covering. Its Jacobian is given by 

d((p,il>) 2u 2u 



d(u,v) ci sin(<£ + tp - a) c \ - [u 2 - c ) ; 

Therefore 



K*f{%) = 4vrr(a) / / f(u) dudv 

Jo JV^l V c i _ ( _ c o) 

= 2r(a) / f{u)— ====== du. 

Jy/^l V C 1 _ ( - C 0) 

To compute A a * A&(/) we set g(x) = * f( x )- Then g is spherical and 

= A a *(A 6 */)(e) 
= (\a*g)(e) 

g(x)d\ a (x) 



S(a) vol. (a), 
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Now 

-2tt /-2tt 



vol. (O b ) = 2r(a) / / dip d^ = 87r 2 r(a). 
Jo Jo 

Therefore 

l*a*IM>(f) = 8n 2 r(a)g(a). 
Substituting from the calculation of (7(a) = A& * /(a) above we obtain 



K * A 6 (/) = (?(a) vol.(a) = 167r 2 r(a)r(6) / /(«) du 

Jsf^l V c i - ( u - c o) 

This completes the proof of part (A). 

Proof of part (C) - Since both / and A a are K-bi- invariant, A a */(x) is i^-bi-invariant 

/ e | A 

and therefore to compute /j, a *f(x) we can assume that x is of the form x = I _t . As 



.0 e - *, 

before let {0 n } be a sequence of spherical functions on G converging weakly to the (singu- 
lar) invariant measure A a on the orbit O a . Applying the polar coordinate decomposition, 
for the convolution A a * f(x) we have 



\**f(z) = / f(v x x )dy 

JOa 

= lim / 6 n (g)f(gx~ l )dg 
n ^°° Jg 



n— >oo 



c lim III n {kia! kz) f (k\a! k2X 1 )8(a!)da! dk\dk2 
' k J k J A 



n^oo 



c lim / / 9 n (a')f(a'kx )S(a')da'dk 
'kJa 



c5(ti) / f(akx x )dk 



K 



Writing M = akx 1 = k\a\k2 , where /ci,/c2 E K = SU(2) and using the coordinates in 

CD 



efj' V-V 71 -P 2 e^ pe"^ 

we compute r in term of ti, ti and fc: 

2 cosh r = Tt(o?) = Tr^a?^ 1 ) = Tr(MM*). (4.3) 
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On the other hand we have 

_ e |(*-*i)^l _ p 2 e W e -|(*i+*) pe -i<P 



M = akx' 1 
Therefore 



Tr( MM*) = 2p 2 cosh(t + h) + 2(1 - p 2 ) cosh(t - i x ) , 
Comparing with (J4.3|) we obtain after a simple calculation 

coshr = coshti cosht — (2p 2 — 1) sinhti sinht 

Now set 

c\ = cosh t\ cosh t, C2 = sinh £1 sinh i 
The function / is spherical so it only depends on the component r and therefore 

„-l> 



\ a *f(x) = c5(t 1 ) / f(akx- L )dk = 05(h) / / / f(r)(2p)dpdipdijj 
Jk Jo Jo Jo 

Now make the change of coordinate coshr = c\ — 02(2 p 2 — 1) and note that, assuming 
that t > t\, r ranges over I t>tl — [t — ti, t + 1\] as p ranges over [0,1]. Substituting in the 
above integral for A a * f(x) we obtain 

2n 2 c f 

A a */(x) = <K*i) / f{r) sinh r dr , (4.4) 

C 2 Jlt, H 

To compute A a * A&(/) we set g(x) = /if, * /(#). Then g is spherical and 

K*W) = A *(A&*/)(e) 
= (A a *#)(e) 

#(x)dA a (x) 



= ^(a) V ol.(0 ), 

Using f)2.2j) one easily obtains 

vo l.(O at ) = c(vol.(5[/(2))) 2 sinh 2 t = 167r 4 sinh 2 t. (4.5) 
Substituting from f)4.5|) we obtain 

A a * A&(/) = 16c7r (sinhti) 2 5f(a). 

Equation (|4.4|) implies 

A a * = 32c 2 7r 6 sinh t x sinh t 2 / /(r) sinh ra?r , 

Jit lt t 2 

which completes the proof of the theorem. ■ 
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5 Discretization and Numerical Results 



The discretization and quantization of conjugacy classes are essentially different issues, 
and the latter is based on the Kirillov orbit method (see [K]). It is of interest in theoretical 
physics and is also relevant to the subject matter of this paper. The dual of the Lie algebra 
of 577(2) is identified with the set of 2 x 2 skew hermitian matrices of trace 0: 

j. / 7 \ / ic a + ib\ 
^ c)= [-a + l b -zcj 

Under conjugation action of G the orbits are the spheres E r = {£(a, b, c) | a 2 +b 2 +c 2 = r 2 }. 
It is customary in physics to assign the symmetric n th power representation p n of £77(2) 
to the sphere E r of area n G Z+. According to the Clebsch- Gordon formula, p n <g) p m 
decomposes as 

n+m 

Pn® Pm — Pk-, where k = m + n mod 2. 

\n—m\ 

On the other hand the Minkowski sum of spheres of radii r\ and r-i in R 3 is precisely the 
spherical shell defined by \r 2 — r±\ < ||£|| < ri + r 2 . The perfect resemblance between sums 
of spheres and the decomposition of the tensor product is carried over to the multiplication 
of conjugacy classes. In fact, we set (note slight change of notation) 

C n = exp(£ji_). 

Let j*- = a mod n and ^ = (3 mod ir. By Theorem 12.11 

C n .C m = exp(E » + Em) 

The representations corresponding to conjugacy classes contained in C n .C m correspond to 
those integers k G n|, m+n] which = m+nmod2. Thus in the range [|m— n|, m+n] 
about half the representation occurring in C n .C m appear in the Clebsch-Gordon formula. 
This is what is meant by the quantization of the product of two conjugacy classes. 

By the discretization of products of conjugacy classes C a and Cp one means the choice 
of iV points on each and the determination of the corresponding empirical measure of 
products of these points. If these iV points are chosen randomly according to the invariant 
measures on the conjugacy classes then the empirical measure of the products converges 
weakly to the density function (J2.1)) by Theorem 12. II and is numerically demonstrated in 
Figures 1 for a typical choice with Np = 2172 points. 

Each conjugacy class C a is naturally equivalent to a copy of £ 2 . It is therefore reason- 
able to investigate the weak convergence of the empirical measure if the points on £ 2 are 
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chosen according to the requirements of Thomson's Problem of distributing points on the 
sphere (see [KS]). This means that the points should be distributed so that the Coulomb 
potential 



is minimized. A variation of this problem for a = 1 was originally posed by J. J. Thomson 
in connection with his investigations of the structure of the atom in 1904. It remains 
unsolved except for a few small values of N, and it has also attracted attention for 
applications to complexity theory [Sin]. The lattice point method makes use of the natural 
embedding of the icosahedron in S 2 and distributes N = I0(mn + m 2 + n 2 ) + 2 points on 
the sphere in such a way that the distribution exhibits a high degree of symmetry and the 
points appear to be "evenly" spaced. It was conjectured in [A] that this distribution will 
provide the solution to Thomson's Problem for a — 1. The polar coordinates method was 
devised in [KS] to achieve the minimum required by Thomson's Problem and numerical 
tests disproved the conjecture in [A] by showing that the (local) minimum achieved by the 
polar coordinates method (where there was symmetry breakdown) was in fact smaller. It 
is therefore natural to test the convergence of the empirical measure of products if the 
discretization is done according the polar coordinates or the lattice point methods. In 
Figures 1 and 2 the convergence of the empirical measure to the density predicted by 
Theorem 12. II is exhibited for Np = 2172 points. The angles </> and ip in the captions refer 
to the conjugacy classes and respectively. In Figures 2 the corresponding empirical 
measures are calculated for the lattice point method and it is noticed that even if the 
measure converges to the density given by (|2.1|) . the convergence is significantly slower. 
Similar conclusion is applicable to the polar coordinates method as shown in Figures 3. 
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<Ht/2 , \|f= 37t/8 , Np=2172 , Method=Polar 
0.45 1 i i 1 i 1 




0.5 1 1.5 2 2.5 3 3.5 Figure 3. 
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